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摘 要 : 为 研究 和 田 河中 下 游 流 域 地 下 水 水 化 学 特征 及 其 演化 规律 ,采用 Piper 三 线 图 、Gibbs 图 、 
PHREEQC 软件 和 数理 统计 等 方法 对 该 地 区 水 化 学 特征 、 主 要 溶质 组 分 来 源 和 演化 规律 开展 了 综 
合 研究 分 析 。 结 果 表 明 :(1) 地 下 水 中 八大 常规 离子 含量 整体 偏 高 ,尤其 是 Cl SOT Na Ca" fü 
HCO; 占 比较 大 。(2) 水 化 学 溶解 组 分 空间 差异 性 明显 , 除 少 部 分 地 区 水 化 学 类 型 表现 为 SO4: Cl- 


Ca* Mg 型 外 ,大 部 分 均 以 SO4*Cl-Na 型 为 主 。 地 下 水 呈 弱 碱 性 ,大 部 分 地 区 水 样 点 水 质 基 本 满足 生 
活 饮 用 水 标准 。 部 分 绿洲 平原 区 受 人 类 生产 活动 影响 ,水 中 NOi 含 量 明显 异常 。(3) xxm TER 


和 阳离子 交换 作用 下 ,地 下 水 中 Nar\Car Mg" CE .SO 等 离子 组 分 主要 物质 来 源 于 岩 盐 方解石、 
白云 石和 石 背 等 矿物 的 溶解 。 在 向 下 游 细 土 平原 和 沙 江 区 径流 过 程 中 , 受 藻 发 浓缩 作用 普遍 控制 
影响 ,水 中 各 离子 浓度 含量 不 断 增 大 。(4) 在 开放 性 较 好 的 系统 环境 中 ,CO 促进 了 名 矿物 的 溶解 
使 得 地 下 水 中 各 离子 含量 不 断 海 解 聚集 ,浓度 增 大 。 进 入 下 游 冲积 平原 区 后 , 受 细 粒 砂 士 介质 层 


生态 环境 保护 提供 理论 依据 。 


阻隔 影响 ,地 下 水 径流 与 交替 强度 减弱 导致 溶 滤 作 用 变 差 ,蒸发 浓缩 作用 逐渐 占 主 导 地 位 ,进一步 
加 大 了 水 中 离子 含量 和 水 化 学 类 型 的 差异 性 。 研 究 结果 可 为 和 田 河 流域 水 资源 合理 
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地 下 水 在 径流 过 程 中 ,不 断 地 与 途经 地 层 发 生 
着 水 - 岩 相 互 作 用 ,进而 改变 着 水 中 的 洲 解 组 分 。 
地 下 水 的 化 学 平衡 状态 和 成 份 演变 主要 受到 地 质 、 
水 文 地 质 条 件 和 气候 等 诸多 因素 的 影响 "1。 

近年 来 ,国内 外 学 者 针对 不 同 区 域 不 同时 间 尺 
度 开 展 了 较为 丰富 的 地 下 水 水 化 学 研究 ,包括 水 化 
学 分 类 ,溶解 组 分 来 源 和 控制 因素 等 不 同 角度 , 研 
究 方法 主要 有 派 珀 (Piper) 三 线 图 "5 、 阿 雇 金 分 类 
法 "Durov 三 线 图 ”等 。 张 勇 军 等 3 在 2018—2021 
年 期 间 , 对 新 疆 喀 什 河 流域 采集 的 324 组 水 样 ( 河 水 
289 组 ,泉水 5 组 , 井 水 30 组 ) 进 行 水 化 学 组 成 及 成 
因 分 析 , 并 探讨 了 流域 水 化 学 的 时 空 变化 特征 和 控 
制 因素 。 杨 锐 等 ”通过 在 新 疆 南 部 和 田 平原 地 区 
进行 地 下 水 设计 取样 240 组 ,利用 水 化 学 测试 结 
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TR ,分 别 在 水 平方 向 和 垂直 方向 上 分 析 了 地 下 水 的 
硬度 分 布 范围 和 空间 影响 因素 ,提出 了 水 - 岩 相 互 
作用 控制 影响 着 地 下 水 硬度 的 分 布 特征 范围 。 李 
玲 等 "在 南 疆 和 田 河 中 下 游 绿 洲 平 原 区 利用 单一 
参数 模型 法 对 地 下 水 质量 评价 及 关键 指标 进行 了 
识别 ,并 分 析 了 该 区 域 地 下 水 质量 在 时 间 与 空间 上 
的 演化 特征 。 除 此 之 外 ,水 文 地 球 化 学 反 向 模拟 越 
来 越 成 为 当前 一 种 使 用 比较 广泛 而 有 效 的 地 下 水 
水 化 学 特征 研究 手段 。 通 过 计算 矿物 相 的 饱和 指 
数 , 模 拟 分 析 地 下 水 与 围 岩 间 的 溶解 .吸附 、 阳 离子 
交换 等 水 岩 相互 作用 过 程 ,能 够 有 效 地 从 定量 的 角 
度 刻画 地 下 水 中 各 种 组 分 的 形成 来 源 、 转 化 方向 和 
变化 趋势 ,在 地 下 水 水 化 学 特征 和 演化 规律 研究 方 
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和 田 河 流域 位 于 新 疆 塔 里 木 盆地 最 南 缘 的 和 
田 县 境内 ,人 口 较 密 集 ,地 下 水 是 主要 的 供水 水 
源 。 前 人 大 量 资料 成 果 表 明 , 和 田 河流 域 水 文 地 质 
工作 程度 较 高 ,目前 地 下 水 水 化 学 研究 成 果 认 识 多 
集中 在 区 域 水 化 学 特征 及 其 影响 因素 分 析 方 面 ,但 
缺乏 对 水 中 化 学 物质 来 源 及 成 因 机 理 的 进一步 研 
究 ,尤其 是 通过 水 文 地 球 化 学 模拟 定量 分 析 描 述 地 
下 水 水 化 学 平衡 状态 和 各 矿物 相 的 相互 转移 量 方 
面 的 研究 仍 属于 空白 区 。 因 此 ,本 文 在 和 田 河流 域 
中 下 游 地 区 21 个 地 下 水 水 样 点 测试 基础 上 ,分 析 影 
响 地 下 水 水 化 学 特征 的 主要 矿物 组 分 以 及 物质 来 
源 , 并 反 向 模拟 地 下 水 与 环境 的 相互 作用 机 制 ,有 
助 于 进一步 深入 了 人 解 地 下 水 的 成 因 和 演化 规律 ,以 
期 为 区 域 地 下 水 资源 合理 开发 利用 和 保护 提供 一 
定 的 理论 依据 。 


1 数据 与 方法 


1.1 研究 区 概况 

人 研究 区 位 于 欧 亚 大 陆 腹地 , 属 典 型 的 内 陆 干旱 
区 ,四 季 多 风沙 ,多 年 蒸发 量 和 降水 量 平均 值 分 别 
为 2480 mm 和 35 mm, 干 旱 指 数 约 71。 地 势 南 高 北 
低 , 整 体 地 貌 单元 依次 表现 为 低 山 丘 陵 区 、 砾 质 平 
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图 1 研究 区 地 下 水 样 点 分 布 及 模拟 路 径 示意 图 


Fig. 1 Distribution diagram of groundwater samples and 


simulated path of the study area 


保 样品 资料 准确 .安全 和 完善 。 水 化 学 溶解 组 分 测 
试 项 目 主要 包括 八大 常规 离子 和 洲 解 性 总 固体 
(TDS) 、 总 硬度 (TH) 等 综合 指标 。 其 中 ,采用 原子 
吸收 光谱 分 析 法 测试 K* 和 Na’ 含 量 , 乙 二 胺 四 乙酸 
(EDTA ) 配 位 滴定 法 测量 Ca* 和 Mg”* 含 量 , 硝 酸 银 容 


原 区 ` 细 土 平 原 区 和 风 积 沙漠 区 ”5 。 第 四 系 地 层 
发 育 良 好 ,成因 类 型 较为 复杂 ,局 部 地 层 厚度 约 在 
千 米 左右 。 地 层 岩 性 结构 由 单一 大 厚度 的 卵 砾 石 、 
砂砾 石 层 逐渐 过 渡 到 砂砾 石 . 中 粗 砂 、 粉 细 砂 等 多 
层 相互 交替 。 地 下 水 类 型 主要 为 第 四 系 松 散 岩 类 
孔 际 潜水 ,从 山 前 冲 洪 积 松散 砂 卵 砾石 平原 区 到 下 
游 的 细 土 平原 区 和 沙漠 绿洲 区 ,地 下 水 位 埋 深 由 大 
变 小 ,变化 范围 约 在 1~30m 之 间 。 地 下 水 主要 补给 
来 源 有 大 气 降水 人 渗 和 绿洲 灌溉 水 人 渗 , 以 及 少量 
的 地 下 水 侧 向 径流 补给 ,最终 多 以 蒸发 蒸腾 、 地 下 
水 侧 向 径流 以 及 农业 开发 等 方式 排泄 "”"。 
1.2 数据 来 源 

依托 中 国 地 质 调查 局 “重点 流域 地 下 水 流 场 调 
查 ” 地 质 调查 项 目 , 存 2021 年 6 一 7 月 丰 水 期 间 , 沿 
和 田 河中 下 游 流域 共 布 设 地 下 水 取样 点 21 个 (图 
1)。 在 野外 使 用 聚 乙 烯 塑料 瓶 采取 水 样 ,过 程 中 严 
格 遵 守 水 质 采样 相关 标准 规范 要 求 。 取 样 时 用 水 
样 冲洗 3 次 后 再 灌 满 并 密封 保存 后 ,及 时 送 至 新 疆 
地 矿 局 第 二 水 文 工 程 地 质 大 队 进 行 检 测 分 析 ,并 确 


量 法 测量 CL 含量 ,硫酸 饥 比 独 法 测量 SOz 含量 , 酸 
碱 滴定 化 学 分 析 法 测量 HCO; 和 CO; 含量, 紫外 光度 
法 测量 NO; 含 量 ,烘箱 直接 烘 干 法 测量 TDS,EDTA 
络 合 滴定 法 测定 TH。 在 样品 测试 过 程 中 ,所 有 水 化 
学 组 分 含量 检 出 限 为 0.05 mg.[ ,并 对 分 析 测 试 结 
果 进 行 阴 阳离子 电荷 平衡 可 靠 性 审查 ,相对 误差 均 
小 于 +5% ,说 明 水 化 学 检测 数据 结果 可 信 。 
13 研究 方法 

在 样品 测试 结果 的 基础 上 ,运用 SPSS 数 理 统计 
软件 对 地 下 水 中 八大 稼 规 离 子 和 部 分 综合 指标 进 
行 描述 统计 以 及 确定 各 指标 间 的 相互 关系 。 利 用 
Piper 三线 图 展现 主要 离子 空间 关系 和 地 下 水 水 化 
学 类 型 。 采 用 Gibbs 图 解 判断 地 下 水 各 化 学 组 成 的 
形成 机 制 和 演化 特征 。 运 用 离子 比例 系数 法 判断 
地 下 水 的 成 因 和 化 学 组 分 的 可 能 矿物 来 源 。 最 后 
通过 PHREEQC 计算 矿物 相 的 饱和 指数 和 反 向 模拟 
路 径 上 的 水 岩 相 互 作 用 过 程 。 另 外 ,采用 Adobe 
Photoshop fll Microsoft Visio 软件 进行 相关 图 件 的 绘 
制 和 处 理 。 


2 结果 与 分 析 


2.1 水 化 学 参数 统计 特征 

人 研究 区 地 下 水 检测 数据 结果 显示 ( 表 1) ,溶解 性 
总 固体 (TDS) 变 化 范围 较 大 ,为 643.5~2295.82 mg: L'', 
平均 值 为 1168.72 mg*L'。 水 中 八大 常规 离子 绝对 
含量 整体 均值 偏 高 ,尤其 是 CI HCO; .Na fI SOF 1H 
比较 大 。 在 质量 浓度 对 比分 析 中 ,阳离子 溶解 组 分 
旺 现 出 Na*>Ca*>Mg”*>K’ 的 一 般 规律 ,平均 浓度 分 
别 为 244.2 mg:L'、80.83 mg .二 、 58.75 mg- L fil 
19.41 mg- L'; 阴离子 则 表现 出 CESHCO:2S0? 2NO; 
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的 特征 ,平均 浓度 分 别 为 352.78 mg-L' 341.5 mg EL, 
235.87 mg .L 404.84 mg- L'', pH f vi FEL YE 7.59- 
8.16 Z [8] ,均值 为 7.88, 整 体 呈 弱 碱 性 。 总 硬度 
(TH) 表 现 出 波动 范围 较 大 ,在 188.4~743.6 mg- L' 
之 间 ,均值 438.87 mg*L'。 根 据 水 的 硬度 表示 方法 
对 研究 区 内 地 下 水 进行 分 类 ,区 内 地 下 水 类 型 主要 
属于 中 硬 水 到 高 硬 水 之 间 "“。 变 异 系数 方面 ,常规 
Bi P K* Na! Ca Me SOF Cr, HCO AK pH, 
TDS TH 值 均 小 于 0.6, 表 明 水 中 离子 浓度 变化 范围 
十 分 稳定 ;而 NO; 离 子 的 变异 系数 较 大 ,为 0.94, 离 子 
含量 最 大 值 为 18.74 mg: L' ,最 小 值 仅 为 0.25 mg L', 
表明 地 下 水 很 有 可 能 受到 人 为 干扰 因素 的 影响 ”。 


R1 地 下 水 水 化 学 参数 统计 


Tab.1 Statistical results of hydrochemical parameters of groundwater 
统计 参数 L' nis ho L' his L' m L' s nun e pH us hos i 
极 差 35.30 428.57 108.77 116.83 638.10 411.52 559.16 18.49 0.57 1652.31 555.20 
最 小 值 9.56 88.70 39.33 21.93 141.80 119.86 147.73 0.25 7.59 643.50 188.40 
最 大 值 44.86 517.27] 148.10 138.76 779.90 531.38 706.89 18.74 8.16 2295.82 743.60 
平均 值 19.41 244.42 80.83 58.75 352.78 235.87 341.50 4.84 7.88 1168.72 438.87 
标准 差 8.32 123.48 33.77 29.61 183.43 92.58 134.30 4.53 0.17 447.83 150.83 
变异 系数 0.43 0.51 0.42 0.50 0.52 0.39 0.39 094 0.02 0.38 0.34 


2.2 水 化 学 类 型 及 空间 分 布 

Piper 三 线 图 常 在 水 文 地 质 学 中 被 广泛 地 用 作 
研究 中 水 样 化 学 成 分 的 主要 离子 组 成 特征 , 它 能 
效 地 将 地 下 水 水 化 学 类 型 进行 分 类 ,也 可 以 用 来 表 
示 水 中 各 成 分 含量 之 间 的 关系 ?59 。 

在 地 下 水 Piper 三 线 图 中 (图 2) ,地 下 水 各 样 点 
集中 分 布 在 2 区 ,水 化 学 类 型 为 SO Cl-Na 型 ,表明 
研究 区 大 部 分 地 下 水 水 样 点 水 化 学 组 分 受到 蒸发 
浓缩 作用 ,这 应 与 图 1 中 大 部 分 水 样 点 位 于 和 田 河 
下 游 的 砾 质 和 细 土 平原 区 有 关 。 其 次 , 少 部 分 水 样 
点 落 在 1 区 ,水 化 学 表现 为 过 渡 的 SO - CI- Ca - Mg 
型 。 仅 有 位 于 和 了 田 河流 域 中 游 砾 石 平 原 区 的 G15 水 
样 点 表现 在 4 区 ,水 化 学 类 型 为 HCO;-Ca.Mg 型 。 
从 各 样 点 水 化 学 阴阳 离子 对 比 上 看 ,左边 阳离子 组 
分 占 比 三 角 图 中 ,水 样 点 多 集中 在 靠近 Na 端 分 布 ， 
说 明 Na 是 主要 优势 离子 。 右 边 阴离子 组 分 占 比 三 
角 图 中 ,水 样 点 多 靠近 在 CL 和 SO, 端 分 布 , 表 明 SOf 
和 CT 含量 占 相 对 优势 。 说 明 研 究 区 内 地 下 水 各 样 
点 中 ,各 自 离子 含量 占 比 存在 明显 差异 性 。 


注 : 表 中 为 各 指标 浓度 的 统计 参数 特征 ;TDS 为 溶解 性 总 固体 ;TH 为 总 硬度 。 下 同 。 
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图 2 地 下 水 Piper 三 线 图 


Fig. 2 Piper trilinear diagram of groundwater 


2.3 离子 来 源 分 析 
地 下 水 水 化 学 溶解 组 分 和 空间 分 布 特征 成 因 
复杂 ,除了 受气 象 水 文 .地 质地 貌 .地 层 岩 性 等 因素 
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控制 影响 外 ,还 往往 与 阳离子 交换 作用 .蒸发 浓缩 
作用 、 水 -岩溶 渡 作 用 、 大 气 降水 作用 和 人 类 活动 影 
响 等 因素 有 关 。 

2.3.1 各 化 学 指标 间 的 相互 关系 采用 主 成 分 分 析 
法 (PCA) 对 水 样 测试 结果 中 Na Ca” Mg” CT SO 
等 11 项 宏 量 组 分 和 化 学 综合 指标 进行 降 维 统计 和 
标准 化 处 理 ( 表 2)。 相 关系 数 和 矩阵 中 数据 显示 LKC 
Na' , Mg" CT SOF 和 TDS 之 间 相 关系 数 都 在 0.7 以 
上 ,表示 地 下 水 中 各 化 学 组 分 指标 间 的 相关 性 良 
好 。 旋 转 特征 值 大 于 1 的 3 个 因子 中 ( 表 3) ,贡献 率 
分 别 占 总 变量 的 55.75% .20.12% 和 9.73% ,三 者 方 
差 累计 贡献 率 达 85.59% ,这 3 个 主 成 分 可 以 反映 出 


总 样本 85.59% 的 信息 量 。 
由 表 4 中 可 以 看 出 ,在 因子 1 中 ,Mg*”*、SO? Na, 


K' .Cl 等 水 化 学 指标 影响 系数 较 高 ,分别 为 0.93、 
0.91 .0.86 .0.86 .0.84, 并 具有 正 向 载荷 量 , 表 明 在 和 
田 河 中 下 游 流域 地 区 内 地 下 水 水 化 学 的 形成 普遍 
受到 蒸发 浓缩 作用 的 影响 控制 。 在 因子 2 中 ,可 以 
明显 看 出 ,Ca” pH 水 化 学 指标 权重 系数 较 大 ,为 主 
影响 因子 成 分 。 反 映 了 地 下 水 在 径流 途经 周围 
岩 土 过 程 中 ,水 岩 涂 滤 作 用 强烈 。 同 时 Ca E f m 
载荷 ,表明 在 细 土 平原 区 内 地 下 水 与 砂砾 土 、 砂 质 
粉 土 等 地 层 发 生 了 强烈 的 阳离子 交换 作用 。 因 子 3 
则 以 NO; 为 主要 影响 成 分 ,结合 表 1 中 NO; 离 子 的 变 
异 系 数 较 大 的 事实 ,说 明 在 下 游 绿洲 平 原 部 分 地 区 
地 下 水 受到 人 类 生产 活动 影响 。 
2.3.2 水 化 学 控制 因素 Gibbs 图 可 用 来 判断 地 下 
水 中 组 分 的 主要 来 源 和 形成 机 制 ,包括 大 气 降 水 、 


表 2 地 下 水 水 化 学 参数 的 相关 系数 徐 阵 


Tab.2 Correlation coefficient matrix of groundwater hydrochemical parameters 


水 化 学 参数 Kt Na* Ca* Mg” CE SOT HCO; NO; pH TDS TH 
K* 1.00 0.72 —0.10 0.80 0.60 0.79 0.75 —0.04 —0.02 0.78 0.58 
Na - 1.00 —-0.22 0.74 0.94 0.82 0.41 0.15 0.13 0.94 0.48 
Ca 一 1.00 0.08 -0.10 0.19 0.36 -0.14 —0.58 0.07 0.63 
Mg” S = S 1.00 0.78 0.75 0.76 —0.07 0.04 0.88 0.81 
cr = = = = 1.00 0.73 035 0.21 0.19 0.93 0.57 
SOF E x E - 一 1.00 0.64 -0.21 -0.17 0.90 0.71 
HCO; - - - - - - 1.00 -0.13 -0.34 0.63 0.78 
NO; - - - - - - - 1.00 0.02 0.06 -0.11 
pH - - - - - - - - 1.00 -0.02 -0.29 
TDS - - - - - - - - - 1.00 0.74 
TH - - - - - - - - - - 1.00 
表 3 相关 和 矩阵 特征 值 及 方差 累计 贡献 率 表 4 因子 载荷 矩阵 
Tab. 3 Eigenvalues of correlation matrix and cumulative Tab.4 Factor loading matrix 
variance contribution rate 因子 
- 水 化 学 参数 
IA 初始 特征 值 1 2 3 
BEEETCTENET TU UUNMETIUEDTT, K 0.86 0.10 -0.15 
1 6.13 55.75 55.75 Na! 0.86 0.43 0.09 
2 221 20.12 75.87 Ca” 0.16 -0.88 0.21 
3 1.07 9.73 85.59 Mg 0.93 0.04 -0.09 
h bus dl 21:09 Cr 0.84 0.40 0.19 
: ~ n mE Sor 0.91 -0.06 -0.14 
6 0.20 1.81 98.86 
7 0.09 0.79 99.65 HCO; 0.77 —-0.41 —0.08 
8 0.03 027 99.91 NO; —0.03 0.34 0.91 
9 0.01 0.07 99.98 pH =0.09 0:76 SM 
10 0.00 0.02 100.00 TDS 0.98 0.15 0.08 
11 0.00 0.00 100.00 TH 0.83 —0.45 0.07 


蒸发 浓缩 和 水 岩 洲 滤 作 用 等 因素 影响 下 的 水 化 学 
演化 特征 "””。 由 图 3 可 知 ,在 图 解 模型 中 右上 角 
第 一 象限 区 域内 ,存在 着 大 部 分 地 下 水 水 样 点 , 绝 
KE Na (Na^ Ca^) fl CE/(CIE--HCO:) EE fet £8 8X 
于 0.6, 最 大 值 可 接近 于 1, 表 明 下 游 靠 近 沙 漠 的 细 
土 平原 绿洲 区 地 下 水 水 化 学 组 分 主要 受到 藻 发 浓 
缩 作用 控制 。 该 结论 与 前 面 主 成 分 分 析 结 果 一 致 ， 
也 符合 研究 区 位 于 气候 干旱 少雨 的 内 陆 干 旱 区 的 
背景 条 件 。 靠 近 中 游 G5、G6、G15 和 G16 水 样 点 基 
本 位 于 水 岩溶 滤 作用 控制 区 域内 ,离子 组 合 比值 一 
般 小 于 0.5 ,说 明 在 该 区 域内 ,水 岩溶 滤 作用 对 地 下 
水 水 化 学 组 分 的 形成 具有 较 大 的 贡献 。 研 究 区 内 
无 任何 水 样 点 位 于 第 二 象限 大 气 降 水 作用 控制 区 
域内 ,说 明 在 和 田 河中 下 游 流 域 地 区 地 下 水 组 分 的 
形成 几乎 与 降雨 没有 任何 直接 关系 。 

2.3.3 离子 比例 及 主要 成 分 来 娠 不 同 成 因 或 不 同 
条 件 下 形成 的 地 下 水 ,各 种 水 化 学 成 分 含量 之 间 的 
比例 系数 往往 有 比较 明显 的 差异 ,因此 可 以 运用 离 
子 比值 法 分 析 地 下 水 的 赋 存 条 件 及 成 因 ,进而 判断 
地 下 水 中 化 学 组 分 的 矿物 来 源 ”。 

CI 是 地 下 水 中 最 稳定 的 离子 ,其 售 量 大 小 一 般 
FITDS 的 含量 呈 线 性 增长 关系 ,因此 常常 可 以 作为 
揭示 地 下 水 水 化 学 演化 过 程 的 示 踪 元 素 。Na* 和 Cl 
两 者 关系 基本 呈 1:1 线 性 增加 (决定 系数 记 =0.98) ,说 
明 地 下 水 在 径流 过 程 中 ,Na* 和 CI 的 来 源 与 岩 盐 (Na- 
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Cl) 的 溶解 密切 相关 (图 4a)。 另 外 , 随 着 Cl 浓度 的 
增加 ,Na .CI 的 摩尔 浓度 比值 YINa)/CD) 仍 保持 在 1: 
1 线 附近 ,进一步 表明 研究 区 内 所 有 地 下 水 水 化 学 
组 分 的 形成 均 受 到 蒸发 浓缩 作用 的 控制 影响 (网 
4b)。y(Ca*+Mg”) 与 YHCO3+SO7 站 离子 组 合 之 间 表 现 
出 明显 线性 关系 , 绝 大 多 数 水 样 点 均 落 在 1:1 平 衡 线 
附近 偏 下 方 (图 4c) ,表明 研究 区 大 范围 内 应 存在 着 
部 分 蒸发 岩 类 (如 石膏 ) 风 化 溶解 参与 地 下 水 中 化 学 
组 分 的 形成 2 2 。 

人 研究 区 内 绝 大 部 分 水 样 y(Na')/y(Ca”) 的 值 均 大 
于 1( 图 4d) ,说 明 区 内 水 化 学 成 分 在 形成 和 演变 过 
程 中 ,普遍 发 生 了 阳离子 交换 吸附 作用 ,这 与 前 面 
主 成 分 分 析 中 结论 一 致 。 碳 酸 盐 的 溶 渡 程度 往往 可 
由 y(SO?+CIT) 和 y(HCO3) 之 间 的 比值 关系 来 确定 *。 
所 有 地 下 水 样 点 均 分 布 在 y(SO?+CI)/y(HC03)=1 比 
值 线 下 方 , 说 明 地 下 水 在 途经 过 程 中 , 均 受 到 碳酸 
岩 或 蒸发 岩 的 深 滤 作用 影响 (图 4e)。 地 下 水 中 SOF 
含量 的 多 少 常常 与 石膏 的 浴 滤 程度 有 关 ,其 浓度 通 
第 可 用 y[(SO7)-(Na'-C1)] 来 衡量 表示 ,而 Ca” 浓 度 则 
往往 由 y(Ca*+Mg”-HCO;) 来 呈现 。 两 者 之 间 的 离子 
组 合 比 值 关 系 可 用 来 判断 地 下 水 SOY 的 物质 成 因 。 
当 两 者 之 比 接近 或 者 达到 1:1 时 ,可 以 基本 判定 地 
层 宕 土 中 的 石膏 矿物 溶解 是 地 下 水 中 SO7 主 要 来 
源 。 绝 大 多 数 地 下 水 水 样 点 分 布 在 1:1 比值 线 附 
近 靠 下 (图 4f) ,表明 区 域内 地 下 水 中 的 S07 含量 
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U vy(Na?) (CE) YINa+Ca) 和 XY(CCL+HCO 分 别 为 其 摩尔 质量 浓度 ;TDS 为 溶解 性 总 固体 。 


图 3 和 田 河 流域 地 下 水 水 化 学 Gibbs 图 
Fig. 3 Gibbs diagram of groundwater hydrochemistry in Hotan River Basin 
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图 4 和 田 河流 域 主 要 离子 及 组 合 关系 图 


Fig.4 Relationships between main ions and their combinations in Hotan River Basin 


小 受 石膏 溶解 程度 的 影响 。 
2.4 水 文 地 球 化 学 模拟 

2.4.1 模拟 路 径 和 “可 能 矿物 相 ” 的 选择 PHREEQC 
反问 地 球 化 学 模拟 对 路 径 的 根本 要 求 是 初始 水 样 
和 最 终 水 样 要 处 于 同一 径流 方向 上 ,在 空间 上 属于 


化 学 溶 渡 组 分 也 是 重要 的 考虑 因素 *”。 利 用 X 射 
线 在 晶体 物质 中 的 衍射 效应 ,对 研究 区 地 下 水 含水 
介质 岩 土 样品 进行 X 衍 射 物 相 分 析 。 结 果 显 示 ,第 
四 系 岩 土 矿物 主要 有 石 富 、 盐 岩 、 文 石 、 方 解 石 \ 日 


上 游 与 下 游 关系 "。 人 研究 区 地 势 南 高 北 低 , 考 虑 地 
形 控 制作 用 影响 ,结合 研究 区 水 文 地 质 条 件 , 本 次 
模拟 在 沿 地 下 水 径流 方向 上 ,从 南 向 北 选取 G5、 
G9, G10 和 G13 4 个 地 下 水 样 点 划分 出 G5 一 G9、 
G9 一 G10 和 G10 一 G13 3 条 代表 性 路 径 (图 1)。 各 取 
样 点 主要 离子 组 分 含量 如 表 5 所 示 。 

在 模拟 过 程 中 ,正确 选取 “可 能 矿物 相 ” 是 水 - 
岩 相互 作用 反应 过 程 中 的 基础 与 关键 ,一 般 是 参考 
研究 区 内 地 下 水 含水 层 介 质 中 的 主要 矿物 成 分 ,水 


云 石 云母 和 高 岭 土 等 。 同 时 ,考虑 到 砾 质 平 原 区 
内 松散 岩 类 孔隙 度 相 对 较 大 ,第 四 系 孔 际 水 位 埋 深 
较 浅 , 水 流 系统 处 于 相对 开放 体系 , 故 在 水 岩溶 解 
反应 的 过 程 中 ,也 应 考虑 C0;, 的 参 入 。 

2.4.2 矿物 的 饱和 指数 和 摩尔 转移 量 饱和 指数 是 
用 来 衡量 水 化 学 组 份 溶解 沉淀 的 重要 指标 ,是 判断 
地 下 水 与 围 岩 矿物 的 溶解 反应 中 处 于 何 种 状态 的 
参数 下 。 通 过 PHREEQC 水 文 地 球 化 学 软件 反 向 模 
拟 ,得 到 可 能 矿物 相 饱 和 指数 和 各 阶段 矿物 转移 量 
计算 结果 ( 表 6、 表 7)。 由 表 6 可 知 ,可 能 矿物 相 中 


*5 路 径 上 样 点 主要 离子 组 分 


Tab.5 Main ion compositions of samples on the simulated path 


编号 K'mg-L'"* Na'/mg:L'' 


Ca"/mg:L"  Mg”/mg L” CI /mg- L” SOi/mg.L' | HCOzmg:L'  TDS/mg.L" pH 

G5 20.48 227.74 76.26 94.70 360.88 231.25 488.36 1259.41 7.87 

G9 23.19 385.31 68.63 66.95 588.47 246.55 256.39 1527.57 7.94 

G10 23.29 302.92 123:22 63.29 410.51 349.96 439.52 1495.24 7.59 
G13 44.86 517.27 102.35 118.56 623.92 531.38 706.89 2295.82 7.65 


白云 石 方解石 和 文 石 的 饱和 指数 值 均 大 于 0 ,说明 
地 下 水 中 这 些 化 学 组 分 正 处 于 过 饱和 状态 ,即将 结 
晶 析 出 。 而 盐 岩 、 硬 石膏 .CO 和 石膏 等 矿物 的 饱和 
指数 均 小 于 0, 表 明 地 下 水 中 这 些 矿 物 组 分 都 处 于 非 
饱和 状态 中 , 仍 将 继续 溶解 。 

结合 表 5 可 知 , 沿 水 样 点 G5 一 G9 一 G10 一 G13 
方向 上 ,地 下 水 中 各 离子 含量 明显 增 大 ,说 明 地 下 


表 6 矿物 饱和 指数 


Tab.6 Mineral saturation index 
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水 在 径流 过 程 中 发 生 了 水 岩 洲 滤 和 蒸发 作用 影 
响 。Na* 和 CI 的 含量 同时 都 是 先 增加 后 接着 减少 最 
后 再 增加 ,而 Ca Mg 与 Na 三 者 含量 基本 呈 相 反 
关系 ,结合 表 7 中 岩 盐 在 各 阶段 上 的 溶解 析出 状态 ， 
说 明 在 模拟 路 径 上 ,伴随 着 明显 的 阳离子 交换 作 
用 ,这 应 与 石膏 高岭土 的 连续 溶解 LE a 方 解 
石和 钙 蒙 脱 石 连续 沉 证 有 关 。 同 时 ,石膏 和 黑 云母 
的 连续 浴 解 分 别 使 得 水 中 SOF- 和 开 ' 离 子 含量 不 断 增 
高 。C0O: 则 发 生 了 先 溶解 后 逸 出 的 过 程 ,说 明 在 开 
放 较 好 的 中 游 冲 洪 积 砂 质 平原 区 中 ,地 下 水 中 C0， 


可 能 矿物 相 饱和 指数 4H 1 — Zu Howe 
fi 土 层 中 的 各 到 容易 被 溶解 。 5 
矿物 名 称 。 分子 式 GS 69 G0 CB fE E a T SED 物 更 容易 被 溶解 bs 
xf CaCO; 0.86 061 070 0.81 JAPA EROA E 5388 JA A BE BA 5] , CO F5) 182 
cum EU 100 075 085 096 指数 明显 呈 增 大 趋势 ,说 明 地 下 水 的 C0, 受到 消耗 ， 
CO; CO; -2.27 -2.61 -2.04 -1.90 _ 
2 2 甘 今 量 氢 3 Ub yE gE c > Er us 
ASE Camco) 245 185 175 232 其 含量 越 来 越 少 , 深 滤 能 力 越 来 越 弱 , 而 蒸发 作用 
石膏 CaS0.2H.0 -1.04 -1.05 -0.67 -0.69 则 不 受 影响 。 在 长 期 不 断 浓缩 作用 下 ,导致 地 下 水 
iA NaCl -5.74 -5.30 -5.57 -5.18 中 各 组 分 浓度 和 TDS 含 量 差异 性 不 断 增 大 。 
表 7 水 文 地 球 化 学 反 向 模拟 结果 
Tab.7 Results of reverse hydrogeochemical simulations 
"m 矿物 转移 量 
B 岩 盐 /mmol*L' A/mm L” 高岭土 /mmol*L'! 钵 蒙 脱 石 /mmol*:L* COxyg ”方解石 nmol SEzifE/mmol*L" 白云 石 /mmol*L" 
G5—G9 «527x10? — «126x10? *8.70x10^ -L05x10*  +1.69x10°  -2.12x10? *7.03x107 -1.76x10^ 
G9—GI0 — -503x10^ . 42.80x10* 42.63x10* -3.7x10* — -1.06x10? — -4.10x10? 42.16x10* -1.64x10? 
Gl0—GI13  4605x10?  -«1.86x10? *6.88x 10^ -8.29x10* . -408x10?  -5.90x10^ +5.56x10* -3.14x10^ 


注 :“+” 值 表示 矿物 发 生 溶解 ;“-" 值 表示 矿物 发 生 沉 演 。 


3 结论 


(1) 和 田 河 中 下 游 流域 地 区 地 下 水 pH 均值 为 
7.88 , 呈 弱 碱 性 ,八大 常规 离子 含量 整体 偏 高 。 阳 离 
子 指标 中 ,浓度 大 小 排序 依次 为 Na’*>Ca”*>Mg”>K’， 
阴离子 则 呈 CI>HCO;>SO0?>N0; 的 关系 。 水 化 学 类 
型 以 SO Cl-Na 为 主 。 在 沿 地 下 水 主要 径流 方 问 
上 ,地 下 水 中 各 离子 浓度 和 TDS 含量 不 断 增 大 。 大 
部 分 地 区 地 下 水 样 点 水 质 总 体 较 好 ,基本 符合 生活 
饮用 水 标准 。 部 分 绿洲 部 分 地 区 受 人 类 生产 活动 
影响 ,NO; 含 量 明显 异常 。 

(2) 地 下 水 在 径流 过 程 中 ,不 断 与 途经 围 岩 发 
生 溶 滤 风 化 作用 和 阳离子 交换 作用 ,从 而 影响 水 化 
学 组 分 的 形成 和 演变 。 下 游 细 土 平原 和 沙漠 区 地 
下 水 水 化 学 特征 和 演化 则 主要 受 蒸发 浓缩 作用 
控制 。 


(3) 模拟 表明 , 岩 土 层 中 岩 盐 .石膏 等 蒸发 岩 的 
溶解 贡献 了 地 下 水 中 Na .CT 和 SOF 的 大 多 数 含 量 ; 
Ca^ Mg 则 是 源 自 于 高 岭 土 和 石膏 等 矿物 的 溶解 ; 
K’ 主 要 源 自 于 黑 云 母 的 连续 溶解 ;C0; 则 是 发 生 了 
先 溶解 后 逸 出 状态 ,说明 C0; 在 开放 性 较 好 的 砂 质 
介质 环境 中 促进 了 各 矿物 的 溶解 ,表现 为 地 下 水 中 
Na .CL 和 SO? 等 组 分 不 断 溶解 聚集 ,离子 浓度 增 大 ; 
但 在 进入 细 土 平原 区 后 ,地 下 水 位 埋 深 变 浅 ,同时 
受到 细 粒 砂 土屋 阻隔 影响 ,地 下 水 径流 速度 放 缓 ， 
水 中 溶解 CO; 含量 减少 ,导致 水 岩溶 滤 作用 变 差 , 薄 
发 浓缩 作用 逐渐 占 主导 地 位 ,使 得 Ca Mg 和 CO 
等 矿物 组 分 多 表现 为 结晶 析出 状态 ,地 下 水 水 化 学 
类 型 逐渐 过 渡 到 SO,:Cl-Na 型 为 主 。 
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Chemical characteristics and evolution of groundwater in the middle 
and lower reaches of Hotan River Basin 
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Abstract: To investigate the chemical characteristics of groundwater and its evolutionary patterns in the middle 
and lower reaches of the Hotan River Basin, Xinjiang, China, this study analyzed 21 groundwater samples from 
the area. It included an examination of the constituents and origins of groundwater solutes and the reverse simula- 
tion of hydrogeochemical processes. The analysis employed Piper trilinear diagrams, Gibbs diagrams, PHREEQC 
software, and mathematical statistics to explore the chemical properties of groundwater, the principal sources of 
solutes, and their evolution in the study region. The results revealed that: (1) High concentrations of eight conven- 
tional ions in the groundwater, with CI, SOF, Na', Ca^, and HCO; being particularly abundant. (2) There was a 
significant spatial variability in the dissolved constituents of the groundwater. The predominant chemical types 
were SO,- Cl-Ca: Mg and SO,: Cl-Na, with the latter being more common across most areas. The groundwater’ s 
alkalinity was generally low, and the water quality in most regions met daily drinking water standards. However, 
in some oasis plains, elevated NO; levels were attributed to human activities. (3) The dissolution of minerals such 
as halite, calcite, dolomite, and gypsum, facilitated by water-rock interactions and cation exchange, was identified 
as the main source of Na', Ca^, Mg”, CI, and SOF ions in the groundwater. During transit to finer soil plains and 
desert areas, ion concentrations increased due to evaporation and concentration processes. (4) In open system con- 
ditions, CO» enhanced the dissolution of various minerals, leading to increased ion concentrations. As groundwa- 
ter flowed into the alluvial plains downstream, the fine sand layer acted as a barrier, reducing the intensity of 
groundwater flow and solute leaching, with evaporation and concentration processes becoming more dominant. 
This study provides a theoretical foundation for the sustainable development and management of water resources 
and environmental protection in the Hotan River Basin. 
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